Observational constraints on hyperons in neutron stars 
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The possibility that neutron stars may contain substantial hyperon populations has important 
implications for neutron-star cooling and, through bulk viscosity, the viability of the r-modes of 
accreting neutron stars as sources of persistent gravitational waves. In conjunction with laboratory 
measurements of hypernuclei, astronomical observations were used by Glendenning and Moszkowski 
[Phys. Rev. Lett. 67, 2414 (1991)] to constrain the properties of hyperonic equations of state within 
the framework of relativistic mean-field theory. We revisit the problem, incorporating recent mea- 
surements of high neutron-star masses and a gravitational redshift. We find that only the stiffest of 
the relativistic hyperonic equations of state commonly used in the literature is compatible with the 
redshift. However, it is possible to construct stifi'er equations of state within the same framework 
which produce the observed redshift while satisfying the experimental constraints on hypernuclei, 
and we do this. The stiffness parameter that most affects the redshift is not the incompressibility 
but rather the hyperon coupling. Nonrelativistic potential-based equations of state with hyperons 
are not constrained by the redshift, primarily due to a smaller stellar radius. 

PACS numbers: 04.40.Dg, 97.60.Jd, 26.60.+C, 95.30.Sf 



I. INTRODUCTION 

The observed masses of neutron stars have been used 
for years to constrain theoretical predictions of the equa- 
tion of state of degenerate matter at high density [3| ■ If 
an equation of state is "soft" (low pressure) at high den- 
sity, the maximum mass of a stable star in general rela- 
tivity is lower than for a "stiff" (high pressure) equation 
of state. The most massive observed neutron star then 
sets a limit on the softness of the equation of state. An 
accumulation of neutron star-cooling observations favors 
(though not decisively) the presence of exotic particles 
such as hyperons in the cores of some neutron stars 
which tends to soften the equation of state. Hyperonic 
couplings in a relativistic mean-field theory of dense mat- 
ter can be constrained to a range of values based on the 
measured properties of A hypernuclei and the maximum 
neutron-star mass 0- This information is in turn useful 
for predictions of the bulk viscosity of hyperonic mat- 
ter P, 1^ , which has important implications for the via- 
bility of r-modes in accreting neutron stars as persistent 
sources of gravitational radiation 0. Spurred by these 
implications and by new observations, we revisit the con- 
straints on hyperonic equations of state. 

First we consider the maximum mass. The most pre- 
cise observations of neutron-star masses come from ra- 
dio pulsars in binaries, which are all measured with 95% 
confidence to be less than 1.5 Mq 0. Accreting neutron 
stars are naturally expected to be more massive, and x- 
ray measurements have long suggested that this is so. 
The best case for decades has been Vela X-1 (4U 0900- 
40), with a most likely mass of about 1.8 M© but with 
1.5 Mq included in the 95% confidence interval .8]. How- 
ever, this measurement is now known to be contaminated 
by oscillations of the high-mass main sequence compan- 
ion ; and while more recent measurements [Tol [lll| can 



claim smaller statistical confidence intervals, they are still 
subject to large and poorly quantified systematic errors. 
Thus 1.5 Mq has remained the constraint for many years. 

However, this is changing. Recent radio observations 
of PSR J0751-hl807 by Nice et al. El IH El yield a 
neutron-star mass greater than 1.6 Mq at the 95% confi- 
dence level. PSR J0751-I-1807 orbits a white dwarf which, 
unlike the main sequence star in Vela X-1, shows no evi- 
dence of oscillations. The orbital period is 6 hours, short 
enough that its decay due to gravitational radiation is 
observable. This provides a post-Keplerian parameter to 
disambiguate the two masses. Marginal detection of the 
Shapiro delay implies intermediate orbital inclination an- 
gles, and disambiguates the inclination angle (somewhat) 
from the mass of the neutron stars (see the figure in 
Ref. E3|)- Also, Ransom et al. jl5] find through measure- 
ments of the periastron advance of the highly eccentric 
orbit that Ter 5 I has 1.68 M© or higher, formally at the 
95% confidence level. Overall this bound is tighter than 
that for PSR J0751-I-1807, but the companion is proba- 
bly a white dwarf and there may be some contamination 
of the relativistic periastron advance by its rotationally 
induced quadrupole. Thus the mass constraint on equa- 
tions of state is now at least 1.6 M© and may be 1.7 M©. 

Another constraint is the measurement of a gravita- 
tional redshift by Cottam, Paerels, and Mendez The 
low-mass x-ray binary EXO0748-676 displayed several 
absorption lines (inferred from multiple x-ray burst spec- 
tra) consistent with a redshift z = 0.35. Estimates of nu- 
merous possible sources of error in the redshift amount to 
a total of no more than 5% E3> implying that equations 
of state should be ruled out if their maximum redshift is 
below about 0.33. 

Qualitatively, it has been stated E3 that softer equa- 
tions of state are disfavored by the redshift; here we make 
that quantitative. Like Glendenning and Moszkowski 
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we constrain the parameters of relativistic mean field 
theory between astronomical observation and hypernu- 
clear experiment. Since there is a fairly large parameter 
space involved, most papers using the results of relativis- 
tic mean-field theory use two canonical parameter sets 
corresponding to a soft equation of state and a stiff equa- 
tion of state. We find that the so-called stiff equation of 
state is actually the softest allowed (marginally) by the 
redshift observation if hyperons are present. However, we 
can and do construct stiffer equations of state that are 
compatible with hypernuclear measurements and consis- 
tent with the redshift observation. The most important 
stiffness parameter as far as the redshift is concerned 
is the hyperon coupling (and thus hyperon population) 
rather than the incompressibility. We also note that 
nonrelativistic potential-based equations of state are not 
greatly constrained by the redshift observation (they can 
all reproduce it due to their smaller stellar radii). The 
new neutron-star masses do not constrain the relativistic 
mean-field equations of state much compared to the red- 
shift, but potential-based models are constrained more 
effectively by the masses than by the redshift. In the 
Appendix we provide tabulations of several relativistic 
mean-field theory equations of state with hyperons that 
satisfy the new observational constraints for a variety of 
nuclear-matter parameters. 



II. EQUATIONS OF STATE 

We consider two types of high-density equation of state 
in this paper. The first is the main focus of the paper, 
and the second is used for comparison to demonstrate 
model dependence. 

The first is based on relativistic mean-field theory and 
is discussed in detail in Ref . m . (We note that other rela- 
tivistic models, such as the relativistic Briickner-Hartree- 
Fock of Ref. iS\ , produce equations of state with qualita- 
tively similar behavior.) Here the low-energy strong nu- 
clear interaction is modeled as the tree-level exchange of 
mesons between baryons (neutrons, protons, and possibly 
hyperons). The starting point is the construction of a rel- 
ativistic Lagrangian, which is a sum of free-particle Dirac 
terms for the baryons and leptons, plus free-particle 
terms for the mesons (scalar a, vector w, and isovec- 
tor p), plus interaction terms including tree-level meson- 
baryon interactions and perturbative self-interactions for 
the cr meson. This makes the theory a phenomenologi- 
cal low-energy effective field theory, although it has the 
advantage of being many-particle and relativistic by con- 
struction so that the sound speed never exceeds the speed 
of light. It also has a small number of parameters which 
can be fit simply to experiment; although this is a mixed 
blessing since the many numbers known from nuclear ex- 
periments must be distilled to a few. Mean-field theory 
also has the disadvantage that it neglects correlations by 
construction. 

Under the assumption that the bulk matter is (on a 



macroscopic scale) static and homogeneous, the fields are 
replaced by their mean values, time and spatial deriva- 
tives vanish, and the Euler-Lagrange equations take a 
form that is relatively simple to solve but is still some- 
what lengthy and thus we do not reproduce it here. It 
is enough to state that the Euler-Lagrange equations in 
this approximation reduce to a set of coupled algebraic 
equations for the lepton and baryon Fermi momenta and 
the meson fields. These are combined with equations for 
generalized /^-equilibrium, electric charge conservation, 
and conservation of baryon number to obtain the Fermi 
momenta and meson fields as functions of, for example, 
the total baryon number density. These are then used to 
construct the pressure and energy density, i.e. the equa- 
tion of state. 

The Euler-Lagrange equations feature five free pa- 
rameters, which under certain assumptions are fit al- 
gebraically to numbers distilled from laboratory mea- 
surements of many finite nuclei: the saturation density, 
binding energy per nucleon and isospin asymmetry coef- 
ficient at saturation density, and the overall incompress- 
ibility K and effective mass m* of nucleons in the nu- 
clear medium. The latter two are difficult to estimate 
from available data and are subject to systematic un- 
certainties, and thus papers using this relativistic set of 
equations of state typically treat a range of values for 
K (240-300 MeV) and m* (0.70-0.80 times the nucleon 
mass m). (The compressibilities are typically lower for 
nonrelativistic models.) 

At roughly twice nuclear density in these equations 
of state, the neutron Fermi momentum is high enough 
to make hyperon production favorable in spite of the 
roughly 200 MeV/c^ mass difference. The hyperons of 
most interest are the A and hyperons, which have 
the lowest masses and therefore are created at the lowest 
densities and occupy the largest fraction of the volume 
of a star. However, the other hyperons S^, S^, and 
'EP also appear in small numbers at the very highest den- 
sities. Hyperons introduce more free parameters. The 
hyperon-meson couplings are assumed to be the same for 
all hyperons but are weaker than the nucleon-meson cou- 
plings by the ratios x^, x^^, and Xp for the three mesons. 
The former two ratios are obtained algebraically from 
the measured binding of A hyperons in nuclear matter 
(double-A hypernuclei) and (more roughly) from hyper- 
nuclear energy levels, resulting in x^ < 0.72 (usually 
taken to be x„ = 0.6) and x^ being determined as a func- 
tion of Xcr and to* . [Note that the first Table in Ref. Q 
contains a typo which is repeated in Ref. 0: The value 
of x^ reading 0.568 should read 0.658, as can be seen by 
solving Eq. (5.59) of the latter reference.] The remaining 
ratio Xp is unconstrained by hypernuclear data, since the 
A is isospin neutral and the relevant E~ hypernuclear 
measurements are highly uncertain. As is standard prac- 
tice, we set Xp = Xo-, although the final equation of state 
is rather insensitive to the precise value 

Numerically, we construct these equations of state us- 
ing the methods of Ref. 1] as functions of the most un- 
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certain parameters m*, K, and Xa- (The low density 
equation of state is the standard BPS model but 
this has little effect on the mass-radius curve which is 
the subject of our work.) For all cases we also use the 
values 0.153 fm""^ for the saturation density, -16.3 MeV 
for the binding energy per baryon at that density, and 
32.5 MeV for the isospin symmetry energy coefficient 
at that density, all as in Ref. iJ. We use seven fidu- 
cial equations of state of this type: Three sets of values 
from Ref. are already in common use in the litera- 
ture: K = 240 MeV and m*/m = 0.78, the softest choice 
which we denote HI; K = 300 MeV and m*/m = 0.78, 
an intermediate choice denoted H2; and K = 300 MeV 
and m* /m — 0.70, the stiffest of these equations of state 
denoted H3. All use Xa — 0.6. We construct the stiffest 
such equation of state compatible with experimental data 
(H4) using K = 300 MeV, m* /ra = 0.70, and x„ = 0.72. 
We construct three others (H5-H7) for extreme values of 
K and to* with x„ just satisfying the astronomical con- 
straints (see below and Fig. |2J). These parameter values 
are summarized in Tabled We also construct relativistic 
mean field equations of state without hyperons by ar- 
tificially setting the hyperon masses to arbitrarily high 
values. These are denoted G1-G7 correspondingly, but 
note that G3 is identical to G4. 

The second type of equation of state is based on de- 
tailed modeling of the potentials observed in laboratory 
nuclei, such as done in Ref. |23|. That paper, denoted 
APR, gives the canonical Schrodinger (non-relativistic) 
model including detailed potentials, two- and three-body 
interactions, and with some relativistic effects in the form 
of perturbations. It has the advantage of using more of 
the known experimental numbers than relativistic mean- 
field theory, including correlations and scattering data. 
However, the fitting to experimental numbers is more in- 
volved, while there are in the end only a few numbers 
that characterize bulk matter and neutron-star struc- 
ture. These equations of state are few-body by con- 
struction and are fundamentally nonrelativistic, result- 
ing in causality violation at high densities (the sound 
speed exceeds the speed of light). Since we expect the 
recent observations to rule out softer equations of state, 



TABLE I: Parameters for seven fiducial hyperonic equations 
of state in relativistic mean field theory. Corresponding equa- 
tions of state without hyperons are denoted G1-G7, but G3 
and G4 are identical. 



Name 


K (MeV) 


m* /m 




HI 


240 


0.78 


0.60 


H2 


300 


0.78 


0.60 


H3 


300 


0.70 


0.60 


H4 


300 


0.70 


0.72 


H5 


300 


0.80 


0.66 


H6 


240 


0.70 


0.67 


H7 


240 


0.80 


0.69 



we also consider a very soft version of this type denoted 
BPAL12 in Ref. The incompressibility of BPAL12 is 
120 MeV, which was known at the time to be much too 
low. It was created explicitly to produce "artificially" 
the softest equation of state compatible with then-known 
neutron star masses (1.45 Mq). Since neither of these 
equations of state includes hyperons, we also consider 
the results of Balberg, Lichtenstadt, and Cook 22], who 
include hyperons in a similar model, and denote their 
equations of state as BLCl (the softer) and BLC2 (the 
stiffcr). 

III. MAXIMUM MASS 

General relativity predicts a maximum mass for a 
star stable to radial perturbations for a given equa- 
tion of state. This is seen by solving the well-known 
Oppcnheimer-Volkoff (OV) equations, which map a curve 
p{p) (pressure as a function of energy density) onto 
a curve M{R) (mass of the star as a function of ra- 
dius). The gravitational mass M generically has a 
maximum, which rules out equations of state that are 
too soft to produce the observed masses. The tight- 
est observational constraint at 95% confidence is now 
1.68 Mq for Ter 5 I though the corresponding 

1.6 Mq for PSR J0751-hl807 [3 may be cleaner. While 
PSR J075H-1807 has a 287 Hz rotation frequency [13, 
this is well below the mass-shedding limit for all equa- 
tions of state and can be shown to increase the OV max- 
imum mass (which assumes no rotation) by no more than 
about 2-3 percent 0|. Ter 5 I rotates at 104 Hz [l5l| . 
and thus its maximum mass is increased by less than 1% 
over the OV value. 

We use a version of the OV equations due to Lind- 
blom 113: 

dm _ ATTp{h)r{hf[r{h) -2m{h)] 

Ih ^ m{h) + 4:Trr{h)3p{h) ' ^ ^ 

dr_ ^ rih)[r{h)-2mih)] 

dh m{h) + ATTr{hfp{h)' ^ ' 

using as independent variable the specific enthalpy 

h{p)= r dp'/[p'+p{p% (3) 

Jo 

(Here G = c = 1.) Unlike the standard OV equations 
for ni(r) (the mass contained within a sphere of radius 
r) and p(r), Lindblom's form does not suffer numerical 
difficulties near the surface of the star, which is simply 
and robustly defined by ft, = 0. We start by picking a 
central enthalpy and evaluating an analytical expansion 
of the equations at a point very close to the center of 
the star (where the equations are singular). We then 
integrate down to ft. = and read off the total mass and 
radius of the star as M = to(0) and R = r(0). 

In Fig. n we plot R{M) for our fiducial relativistic 
mean-field equations of state, with hyperons (bottom) 
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and without hyperons (top). All equations of state with- 
out hyperons are consistent with 1.68 Mq. The softest 
one with hyperons (HI) has a maximum mass of 1.55 Mp 
and is ruled out by the new pulsar observations 14, 1511. 
H2 is nominally inconsistent with the 95% confidence 
limit of 1.68 A/© for Ter 5 I. However, in practice H2 
cannot be ruled out by this observation and must be 
considered marginally consistent because rotation (not 
included in the OV model) can account for most of the 
0.02 Mq difference and changing the confidence level very 
slightly from 95% would bring it within the limit. At the 
68% confidence level for PSR J0751+1807 (1.8 M©), 
only the stiffest equation of state with hyperons (H3) 
is marginally allowed, while again all equations of state 
without hyperons are allowed. 

Of the potential-based equations of state (plotted 
in 0j 22]), BPAL12 is firmly ruled out with a maximum 
mass of 1.45 Mq (at a radius of 9.0 km). This is not too 
surprising, since BPAL12 was deliberately constructed 
with an artificially low incompressibility K — \2Q MeV 
as an extreme example. The extremely soft example for 
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FIG. 1: Oppenheimer-VolkoflF mass-radius curves for the fidu- 
cial relativistic mean field equations of state with hyperons 
(bottom) and without hyperons (top). The vertical line is 
the observational 95% confidence limit 1.68 Mq from Ter 5 I. 
HI is ruled out by the observation. H2 is marginally consis- 
tent, allowing for a small increase in maximum mass due to 
rotation and the imprecision of the constraint. 



this type of equation of state should now be BPAL21, 
which has a maximum mass of 1.67 Mq for a nonrotat- 
ing star (the 95% confidence limit from Ter 5 I) at a 
radius of 9.2 km. The hyperonic BLCl equation of state 
has a maximum mass of 1.55 Mq, which is also firmly 
ruled out. The stiffer BLC2 has a maximum mass of 
1.75 M©, which is compatible with the 95% confidence 
limit of Ter 5 I. APR stars have a maximum mass of 
2.2 Mq, compatible with all constraints. 



IV. GRAVITATIONAL REDSHIFT 

General relativity also predicts a redshift for photons 
leaving the surface of a star with a strong gravitational 
field. For a nonrotating star the redshift z obeys the 
relation 



2GM\ 



-1/2 



(4) 



Since R decreases with M as M approaches its maxi- 
mum for a stable star, z generically has a maximum for 
the maximum-mass star and an observation can rule out 
equations of state which cannot produce a strong enough 
redshift. Cottam, Paerels, and Mendez jl6j| have such an 
observation, a gravitational redshift z = 0.35 obtained 
by identifying several absorption lines in spectra con- 
structed from multiple bursts from the low-mass x-ray 
binary EXO0748-676. Due to the number of consistent 
lines the result is robust, although there may be errors 
at the few percent level Pulsations from a more 

recent x-ray burst have inferred a rotation frequency of 
45 Hz for the neutron star |2^. At this frequency rota- 
tional corrections to the redshift should be a fraction of a 
percent fl^ and the nonrotating approximation suffices. 

The redshift as a function of mass is plotted for 
the fiducial relativistic mean-field equations of state in 
Fig. Eland compared to the observational constraint from 
EXO0748-676. Without hyperons, all equations of state 
are consistent with z — 0.35 for masses greater than the 
1.4 Mq typical of previous measurements, consistent with 
the suspected higher masses of accreting stars in x-ray bi- 
naries. With hyperons, only the stiffest of the usual equa- 
tions of state (H3) is marginally consistent with z = 0.35. 
Therefore we favor using H3 as the new "soft" equation of 
state of this type and H4 as the stiffest. In fact, varying 
values of the nuclear incompressibility K and nucleon 
effective mass m* allow for several "soft" equations of 
state marginally consistent with the redshift, as shown 
in Fig.O Although H3 nominally has z < 0.34, it should 
be considered marginally consistent because there may 
be measurement errors of order 5% and, more im- 
portantly, the maximum redshift is extremely sensitive 
to the equation of state at several times nuclear density. 
For example, artificially excluding all hyperons but the A 
and S~ raises the maximum redshift by 0.03 (a 10% cor- 
rection), even though the populations of those hyperons 
are very small. 
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FIG. 2: Gravitational redshift vs. mass for relativistic mean 
field theory equations of state with hyperons (bottom) and 
without hyperons (top). The horizontal line is 2 = 0.35 mea- 
sured for EXO0748-676. HI and H2 are ruled out, even al- 
lowing for the estimated measurement errors. H3, formerly 
considered the "stiff" equation of state, is actually the softest 
compatible with the redshift. H5 through H7 barely satisfy 
the redshift constraint by construction (see the next Figure). 



Of the potential-based equations of state, BPAL12 is 
marginal with z < 0.36 while aU the others (including 
those with hyperons) easily meet the observational con- 
straint, even when the maximum masses are similar to 
the relativistic mean-field theory models. This is be- 
cause the maximum- mass stars have 11 km radii in the 
relativistic models and 9 km radii in the nonrelativis- 
tic ones. Radii in general are approximately determined 
by the pressure near nuclear density |2^ . which in rel- 
ativistic mean-field theory is about twice what it is for 
potential-based models. Physically this has the simple 
explanation that most of the matter in the neutron star is 
within a factor of two of nuclear density, and so the pres- 
sure at higher densities matters less for the typical radius 
(though it is important for the maximum mass). The fact 
that the redshift does not constrain the potential-model 
equations of state suggests that their low pressure near 
nuclear density may be favored (in the sense that there 
is more unconstrained parameter space). However, at 
high densities these models violate causality, which then 
favors the relativistic mean field models at high density. 




K (MeV) 300 m*/m 



FIG. 3: Relativistic mean field theory equations of state 
can be described by incompressibility K (in MeV), effective 
nucleon mass m* /m, and scalar meson-hyperon coupling Xa 
(see text). The dark surface marks those equations of state 
with a maximum redshift of 0.35. Equations of state below 
the surface are incompatible with the observed redshift of 
EXO0748-676. While H3 is below the surface, it is within the 
estimated 5% error bar and should be considered marginally 
allowed. Equations of state above the lighter surface are com- 
patible with the 95% confidence limit on the mass of Ter 5 I. 
The points corresponding to our fiducial equations of state 
are indicated. 



Rather than using a set of fiducial equations of state, 
one can invert the problem to ask "Given an observation 
of z = 0.35, what parameters in relativistic mean-field 
theory with hyperons are compatible with it?" As stated 
in Sec. II, the three main parameters are Xu, K, and 
m* /m. If we put them all at the stiff end of their al- 
lowed ranges consistent with nuclear and hypernuclear 
experiment, we obtain equation of state 114 with a max- 
imum redshift z = 0.43. If we soften the equation of 
state by lowering Xa- (and thus x^ and Xp) while keep- 
ing the other parameters fixed at the 114 values, we get 
z = 0.35 at x^ — 0.61. Similarly, softening the equa- 
tion of state by lowering K or raising m*/m results in 
bounds of K > 210 MeV or m*/m < 0.84 respectively, 
both of which are less stringent than the experimental 
bounds 0. A plot showing the redshift constraint sur- 
face (and Ter 5 I mass constraint surface) in the three- 
parameter space is shown in Fig|3| The boundaries of the 
cube correspond to the parameter ranges inferred from 
experiment, except for Xa where the experimental lower 
bound is far below the redshift constraint surface. The 
redshift seems to be most sensitive to the hyperonic cou- 
pling parameter x^, and fairly insensitive to the incom- 
pressibility K which is traditionally considered the mea- 
sure of stiffness. This can be seen from the fact that the 
redshift constraint surface in Fig. O is fairly flat and that 
Xa is the only difference between H3 and 114. (The mass 
constraint surface is tilted, i.e. more correlated with K 
and TO*.) The physical explanation is that the main hy- 
peron interaction at high densities is the repulsion (rep- 
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resented in this framework by the vector meson w), and 
thus increasing the couphng constants decreases the hy- 
peron population of a given star. Hyperons, and any 
other new degrees of freedom, soften the equation of state 
at high densities and thus reduce the maximum mass and 
redshift. 

There is also a very recent discussion of a measured 
z = 0.4 in the x-ray binary 4U 1700+24 by Tiengo et 
al. but it is very tentative. This redshift comes from 
one spectral line, which is probably better explained by 
z = 0.012 (implying emission well away from the surface 
of the star) , and there are no other spectral features con- 
sistent with z = 0.4. H4 would still be compatible with 
such a redshift, but would be fairly marginal. 

V. CONCLUSION 

We have compared equations of state for hyperon stars 
with new astronomical observations of mass and gravi- 
tational redshift. Nonrelativistic potential-based models 
are not greatly constrained by the new observations. Rel- 
ativistic mean-field theory models, however, are tightly 
constrained by the observed gravitational redshift. In 
fact, the stiffest of these models commonly used in the 
literature (which we denote H3) is so soft as to be only 
marginally compatible with the observation. When the 
full range of parameters consistent with experiments on 
hypernuclei is considered, there are still many such equa- 
tions of state allowed and the hyperon coupling parame- 
ter is found to be the main one determining the redshift. 
As a consequence we advocate that future studies in- 
volving these models use H3 and a new set of canonical 
parameter values which we denote as equations of state 
H4-H7 (see the Appendix). If through further observa- 
tions the 95% confidence limit on the mass of the neutron 
star in PSR J0751-I-1807 is narrowed to the present 68% 
confidence limit, it would rule out all but H4. 

Moving away from the details of a specific model, the 
general physical result is this: The presence of hyperons 
in neutron stars is constrained but not ruled out by the 
gravitational redshift observation (and to a lesser extent 
by high mass observations). In general the equations of 
state that survive are stiffer than the range previously 
considered in the literature. This means that the hyper- 
ons are less numerous, reducing for example the effect 
of enhanced cooling and bulk viscosity, which is a sub- 
ject for future work It may also be useful to con- 
sider phenomenological equations of state that behave 
like potential-based models near nuclear density but like 
relativistic mean-field models at higher densities. 
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APPENDIX 

For many purposes the detailed microscopic proper- 
ties of a matter model arc unnecessary, and all that is 
desired is a tabulation of pressure, energy density, and 
baryon number density. In Tables II, III, IV, V, and VI 
we provide these for equations of state H3~H7. The low 
density {n < 0.3) parts of some of these (before hyper- 
ons or other strange matter appear) are the same as in 
Ref. . We duplicated the procedure of Glendenning H 
from the beginning using constants from Ref. j30| , which 
leads to some discrepancies in the third or fourth signif- 
icant figure. 

All equations of state but H6 are given up to a baryon 
density n = 1.2 fm~'^, which is more than sufficient for 
stable nonrotating stars. We stop H6 at n = 0.81 be- 
cause at high densities the effective mass of the proton 
becomes negative. This indicates a limitation of the La- 
grangian, which was posited as a low-energy effective the- 
ory. In practice this is not an issue since n = 0.81 is 
almost the central density of the maximum-mass non- 
rotating H6 star. We find that extrapolating H6 beyond 
n = 0.81 under a wide range of assumptions only changes 



TABLE II: Baryon number density n (fm ^), energy density 
p (erg/cm'^), and pressure p (dyn/cm^) for H3. 





n 




P 




P 




n 




P 




P 





03 


5 


041e-fl3 


3 


581e-F31 





63 


1 


221e-hl5 


1 


842e-h35 





06 


1 


009e-fl4 


3 


148e-F32 





66 


1 


290e-hl5 


2 


024e-h35 





09 


1 


518e-H4 


1 


237e-F33 





69 


1 


359e-M5 


2 


217e-H35 





12 


2 


031e+14 


3 


124e-h33 





72 


1 


430e-Hl5 


2 


407e-H35 





15 


2 


549eH-14 


6 


126eH-33 





75 


1 


501e-H5 


2 


597e-f35 





18 


3 


073e+14 


1 


040eH-34 





78 


1 


573e-H5 


2 


791e-H35 





21 


3 


604e-Hl4 


1 


624e+34 





81 


1 


646e+15 


2 


991e-H35 





24 


4 


144eH-14 


2 


389e-H34 





84 


1 


720e+15 


3 


197e-f35 





27 


4 


694e+14 


3 


360e-H34 





87 


1 


794e-fl5 


3 


408e-f35 





30 


5 


263e-hl4 


4 


261e-h34 





90 


1 


870e-Hl5 


3 


596e-f35 





33 


5 


846e-Hl4 


5 


122e-h34 





93 


1 


946e+15 


3 


773e-f35 





36 


6 


440e-fl4 


6 


107e-f34 





96 


2 


023e+15 


3 


950e-f35 





39 


7 


044e-Hl4 


7 


236e-f-34 





99 


2 


lOOe+15 


4 


125e-H35 





42 


7 


658e-H4 


8 


417e-H34 


1 


02 


2 


178e-H5 


4 


301e-H35 





45 


8 


282e-Hl4 


9 


588e-H34 


1 


05 


2 


256e+15 


4 


478e-H35 





48 


8 


914e-Hl4 


1 


083e-f-35 


1 


08 


2 


336e+15 


4 


657e-f35 





51 


9 


556e-H4 


1 


215e-h35 


1 


11 


2 


416e-Hl5 


4 


837e-H35 





54 


1 


021e-Hl5 


1 


357e-H35 


1 


14 


2 


496e-fl5 


5 


019e-f35 





57 


1 


087e-H5 


1 


508e-H35 


1 


17 


2 


577e+15 


5 


202e+35 





60 


1 


153e-fl5 


1 


670e-f35 


1 


20 


2 


659e-fl5 


5 


386e-f35 



7 



TABLE III: Same as the previous Table, but for H4. 



the maximum mass of a stable star by 1%. 



n 



P 



0.03 
0.06 
0.09 
0.12 
0.15 
0.18 
0.21 
0.24 
0.27 
0.30 
0.33 
0.36 
0.39 
0.42 
0.45 
0.48 
0.51 
0.54 
0.57 
0.60 



5.0416+13 
1.009e+14 
1.518e+14 
2.031e+14 
2.549e+14 
3.073e+14 
3.604e+14 
4.1446+14 
4.694e+14 
5.256e+14 
5.839e+14 
6.438C+14 
7.0506+14 
7.676e+14 
8.3156+14 
8.967e+14 
9.631e+14 
1.031e+15 
1.1006+15 
1.1706+15 



3.5816+31 
3.1486+32 
1.2376+33 
3.124e+33 
6.1266+33 
1.0406+34 
1.624e+34 
2.3896+34 
3.3606+34 
4.5616+34 
5.711e+34 
6.924C+34 
8.3136+34 
9.883e+34 
1.1496+35 
1.3216+35 
1.508e+35 
1.710e+35 
1.9286+35 
2.1646+35 



n 



P 



0.63 
0.66 
0.69 
0.72 
0.75 
0.78 
0.81 
0.84 
0.87 
0.90 
0.93 
0.96 
0.99 
1.02 
1.05 
1.08 
1.11 
1.14 
1.17 
1.20 



1.2426+15 
1.3156+15 
1.390e+15 
1.4666+15 
1.543C+15 
1.6216+15 
1.7016+15 
1.782C+15 
1.8656+15 
1.949e+15 
2.034e+15 
2.1206+15 
2.2086+15 
2.2966+15 
2.3866+15 
2.4766+15 
2.569e+15 
2.662e+15 
2.7566+15 
2.8526+15 



2.4166+35 
2.6856+35 
2.9726+35 
3.2606+35 
3.5566+35 
3.8636+35 
4.1826+35 
4.512e+35 
4.8566+35 
5.1916+35 
5.517e+35 
5.8476+35 
6.1836+35 
6.5256+35 
6.875C+35 
7.2326+35 
7.597e+35 
7.969e+35 
8.3466+35 
8.7276+35 
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TABLE IV: Same as the previous Table, but for H5. 



n 




P 




P 




n 




P 




P 


0.03 


5 


045e+13 


4 


568C+31 





63 


1 


207e+15 


1 


720e+35 


0.06 


1 


OlOe+14 


2 


739e+32 





66 


1 


274e+15 


1 


879e+35 


0.09 


1 


518e+14 


1 


123e+33 





69 


1 


343e+15 


2 


045e+35 


0.12 


2 


032e+14 


2 


929e+33 





72 


1 


412e+15 


2 


219e+35 


0.15 


2 


549e+14 


5 


788e+33 





75 


1 


482C+15 


2 


401e+35 


0.18 


3 


073e+14 


9 


754e+33 





78 


1 


552e+15 


2 


591e+35 


0.21 


3 


602e+14 


1 


494e+34 





81 


1 


624e+15 


2 


789e+35 


0.24 


4 


140e+14 


2 


139e+34 





84 


1 


697C+15 


2 


995e+35 


0.27 


4 


686e+14 


2 


916e+34 





87 


1 


770e+15 


3 


209e+35 


0.30 


5 


241e+14 


3 


827e+34 





90 


1 


844e+15 


3 


431e+35 


0.33 


5 


806e+14 


4 


874e+34 





93 


1 


919e+15 


3 


660e+35 


0.36 


6 


388C+14 


5 


840e+34 





96 


1 


995C+15 


3 


898e+35 


0.39 


6 


982e+14 


6 


819e+34 





99 


2 


071e+15 


4 


136e+35 


0.42 


7 


586e+14 


7 


883e+34 


1 


02 


2 


149e+15 


4 


377e+35 


0.45 


8 


199e+14 


9 


040e+34 


1 


05 


2 


227c+15 


4 


623C+35 


0.48 


8 


822e+14 


1 


029e+35 


1 


08 


2 


306e+15 


4 


875e+35 


0.51 


9 


454e+14 


1 


156e+35 


1 


11 


2 


386e+15 


5 


132e+35 


0.54 


1 


OlOe+15 


1 


287e+35 


1 


14 


2 


466e+15 


5 


396e+35 


0.57 


1 


074e+15 


1 


424e+35 


1 


17 


2 


548e+15 


5 


666e+35 


0.60 


1 


140e+15 


1 


558e+35 


1 


20 


2 


630e+15 


5 


942e+35 



TABLE V: Same as the previous Table, but for H6. 



n 




P 




P 




n 




P 




P 


0.03 


5 


038e+13 


2 


458e+31 





45 


8 


220e+14 


8 


987e+34 


0.06 


1 


009C+14 


3 


597C+32 





48 


8 


846e+14 


1 


026e+35 


0.09 


1 


517e+14 


1 


380e+33 





51 


9 


481e+14 


1 


166e+35 


0.12 


2 


031e+14 


3 


308e+33 





54 


1 


013e+15 


1 


320e+35 


0.15 


2 


549e+14 


6 


183e+33 





57 


1 


078e+15 


1 


489e+35 


0.18 


3 


073e+14 


1 


OOle+34 





60 


1 


145e+15 


1 


674e+35 


0.21 


3 


603e+14 


1 


513e+34 





63 


1 


212e+15 


1 


875e+35 


0.24 


4 


139e+14 


2 


160e+34 





66 


1 


281e+15 


2 


090e+35 


0.27 


4 


685C+14 


2 


974C+34 





69 


1 


350e+15 


2 


299e+35 


0.30 


5 


241e+14 


3 


933e+34 





72 


1 


421e+15 


2 


512e+35 


0.33 


5 


816e+14 


4 


709e+34 





75 


1 


493e+15 


2 


732e+35 


0.36 


6 


402e+14 


5 


588e+34 





78 


1 


565e+15 


2 


959e+35 


0.39 


6 


998e+14 


6 


621e+34 





81 


1 


639e+15 


3 


188e+35 


0.42 


7 


604e+14 


7 


807e+34 















TABLE VI: Same as the previous Table, but for H7. 



n 




P 




P 




n 




P 




P 


0.03 


5 


039e+13 


3 


353e+31 





63 


1 


195e+15 


1 


621e+35 


0.06 


1 


009e+14 


3 


585e+32 





66 


1 


261e+15 


1 


777e+35 


0.09 


1 


518e+14 


1 


329e+33 





69 


1 


328e+15 


1 


942e+35 


0.12 


2 


031e+14 


3 


151e+33 





72 


1 


396e+15 


2 


116e+35 


0.15 


2 


549C+14 


5 


842e+33 





75 


1 


465G+15 


2 


298e+35 


0.18 


3 


072e+14 


9 


426e+33 





78 


1 


535e+15 


2 


488e+35 


0.21 


3 


601e+14 


1 


401e+34 





81 


1 


605e+15 


2 


715e+35 


0.24 


4 


136e+14 


1 


966e+34 





84 


1 


677e+15 


2 


897e+35 


0.27 


4 


678e+14 


2 


645e+34 





87 


1 


749e+15 


3 


114e+35 


0.30 


5 


228e+14 


3 


442e+34 





90 


1 


822e+15 


3 


339e+35 


0.33 


5 


787e+14 


4 


362e+34 





93 


1 


896e+15 


3 


572e+35 


0.36 


6 


356C+14 


5 


349C+34 





96 


1 


971C+15 


3 


812e+35 


0.39 


6 


942e+14 


6 


252e+34 





99 


2 


047e+15 


4 


059e+35 


0.42 


7 


537e+14 


7 


227e+34 


1 


02 


2 


124e+15 


4 


309e+35 


0.45 


8 


140e+14 


8 


292e+34 


1 


05 


2 


201C+15 


4 


564C+35 


0.48 


8 


752e+14 


9 


454e+34 


1 


08 


2 


279e+15 


4 


825e+35 


0.51 


9 


374e+14 


1 


071e+35 


1 


11 


2 


358e+15 


5 


093e+35 


0.54 


1 


OOOe+15 


1 


199e+35 


1 


14 


2 


438e+15 


5 


369e+35 


0.57 


1 


064e+15 


1 


332e+35 


1 


17 


2 


519e+15 


5 


651e+35 


0.60 


1 


129e+15 


1 


473e+35 


1 


20 


2 


601e+15 


5 


941e+35 



